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Abstract 

Microwave emissivity models of sea ice are poorly validated empiri- 
cally. Typical validation studies involve using averaged or stereotyped 
profiles of ice parameters against averaged radiance measurements. 
Measurement sites are rarely matched and even less often point-by- 
point. Because of saline content, complex permittivity of sea ice is 
highly variable and difficult to predict. Therefore, to check the valid- 
ity of a typical, plane-parallel, radiative-transfer-based ice emissivity 
model, we apply it to fresh water ice instead of salt-water ice. Radiance 
simulations for lake ice are compared with measurements over Lake 
Superior from the Advanced Microwave Scanning Radiometer on EOS 
(AMSR-E). AMSR-E measurements are also collected over Antarctic 
icepack. For each pixel, a thermodynamic model is driven by four years 
of European Center for Medium Range Weather Forecasts (ECMWF) 
reanalysis data and the resulting temperature profiles used to drive 
the emissivity model. The results suggest that the relatively simple 
emissivity model is a good fit to the data. Both cases, however, show 
large discrepencies whose most likely explanation is scattering both 
within the ice sheet as well as by cloudy atmospheres. Scattering is 
neglected by the model. Further work is needed to refine the scattering 
component of ice emissivity models and to generate accurate estimates 
of complex permittivities within sea ice. 
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1 Introduction 



In Mills and Hevgster ( 2011bl ). a plane-parallel radiative transfer (RT) model 
was used to simulate emitted microwave brightness temperatures at L-band 
(1.4 GHz) microwave frequencies. The question addressed in this paper is 
simply the following: is such a model a reasonable approximation to reality? 
Can factors not accounted for in the m odel, such as ice ridging an d surface 
scattering, be neg lected? The results iMills and Hevesterl ^201 lbl ) suggest 
that while the effect of ice ridging is small, it is not insignificant. 

Because of the simplicity of the model and because effective permittiv- 
ities, required as input to the model, are notoriously diffic ult to estimate 
for saline ice (|Ulabv et all I198& IMills and Hevgsterl . 1201 lbl ) , such a model 
would be a good candidate for inverse methods that aim to retrieve physical 
properties of sea ice beyond mere concentration. The idea is to retrieve those 
quantities most relevant to the model, namely the complex permittivities, 
rather than going one step further back and trying to retrieve ice bulk and 
microstructural properties. The model is only appropr iate to lower (below 
25 GHz) frequencies where volume scattering is wea k ( Mills and Heygsterl 



201 lbl : barber et all Il998l : IVant et all Il978l : Istogrvnl . Il98fil : I.Tohnseil 



The functional dep endence of permitt ivity with frequency in this range also 

19781 ) . thus it may be possible to assume 



tends to be weak (Vantetal 



constant permittivity across all frequencies. 

With the launch of the new SMOS instrument operating at L-band, the 
possibilities for such a retrieval expand considerably. With SMOS, there is 
some potential to retrieve ice thickness. In the event that ice is too opaque 
for this, however, another good candidate for retrieval is snow-thickness. 
Dry snow is nearly transparent at L-band but becomes progressively more 
opaque at higher frequencies (see Equation 0]) . 

The motivation for this study is two-fold: first, there is a near absence 
of direct, point-by-point validation studies of ice emissivity models of this 
type, that is, by taking ice cores, feeding the measured profiles into an emis- 
sivity model and comparing the results to radiance measurements taken at 
corresponding locations. Most validation studies consist of comparing aver- 
aged measured brightness temperatures with simulations based on averaged 
or idealized vertical ice profiles. Second, the complex permittivity of pure 
ice is much less variable than that of saline ice, making it easier to predict. 
Because of the very low value of the imaginary part, which determines the 
attenuation, the effect of ice thickness on the signal over thin ice, i.e., over 
lake ice, will be relatively small. 

The procedure will be to simulate microwave emission over freshwater 
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ice- Lake Superior and the Antarctic ice cap-and compare these with satel- 
lite measurements from the Advanced Microwave Scanning Radiometer on 
EOS (AMSR-E). 



2 Models 

2.1 Emissivity model 



The rad iative transfer model used here is described thoroughly in lMills and Hevgster 



(j2011bh so only the briefest treatment will be given here. A plane parallel ge- 



ometry is assumed making it possible to solve for upwelling and downwelling 
radiances along a single line-of-sight: 

Ti t -7i(l - Ri)T i+1 t -nRiT, | = (1 - n)T, (1) 
T % | - n (l - Ri-x)^ 1 -nRi-xTi f = (l - n )T t , (2) 

where Tj is the physical temperature of the ith layer, Tj f and T, 4- are the 
upwelling and downwelling radiances (as brightness temperatures) respec- 
tively. If there are discontinuous interfaces between the layers, the reflection 
coefficients, Ri can be calculated from the Fresnel equations. The transmis- 
sion coefficients, Tj, is derived: 

/ aiAzA 

n = exp — , (3) 

V cos 9i J 

from the layer thickness, Az, the transmission angle, 0{ (calculated from 
Snell's law), and the attenuation coefficient, af. 

^ v ■ ( A\ 

OLi = imagnj, (4) 

c 

which, in turn, depends on the frequency, v and the imaginary part of the 
refractive index, n« = ^/Ti. c is the speed of light and ej is the relative 
permittivity. 

To run the model, we will need relative permittivities as input. The real 
part for pure ice tends to be fairly const ant at aro u nd 3. 15. The following, 
temperature-dependent model is given in iMaetzlerl (j2006l ): 

e' pi = 3.1884 + 9.1 x 1(T 4 T (5) 

where T is the physical temperature in degrees Celsius. Ice is almost a per- 
fect dielectric: the imaginary permittivity at microwave frequencies is very 
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small, generally less than 0.01. We use the model from iHufford (|l99lh to 
model the imaginary permittivity of pure ice as a function of temperature 
and frequency. These estimates will suffice for lake ice, however the density 
of ice pack in the Antarctic varies between less than 500 kg/m 3 at the surface 
to over 900 kg/m 3 or close to the density of pure ice, at tens of metres depth. 
Since the ice pack is typically granular, we use a mi xture model for spher- 
ical i nclusions to calculate the effective permittivity (jSihvola and au Kongl . 
1988h : 

, 3/(e 2 -ei)ei/(ei + 2e 2 ) 
6 =€l+ l-f(e 2 -e 1 )/(e 2 + e 1 ) (6) 
where e± is the complex permittivity of pure ice, e 2 = 1 is the permittivity 
of air and / = 1 — p/p p i is the relative volume of air, which we can calculate 
from the density of the snowpack, p, versus the d ensity of pure ice, 

We use the following exponential model from iRist et al.l ( 2002 ) for the 
density of snowpack cLS EL function of depth: 



where z is depth. 



■ = 918. - 539.exp(-z/32.5) 



(7) 



2.2 Thermodynamic model 

For the 6GHz channel, the most important determinant of the microwave 
signature of Antarctic ice pack is the temperature and our results will show 
this. For higher frequencies, the temperature becomes less important, but is 
still nonetheless significant. Because of the high penetration depth of low- 
frequency microwaves in pure ice, temperatures quite deep in the ice will 
affect the signal. They will also vary quite significantly from those on the 
surface because of thermal conduction lag. Thermal conduction is modelled 
in one dimension with the diffusion equation: 

at pc p dz 2 { ' 

where k is the thermal conductivity of ice, C p is its heat capacity and po 
its maximum density. The top layer is forced from European Centre for 
Medium Range Forecasts (ECMWF) fields of surface air temperature, T a , 
surface wind speed, v, humidity, q, and cloud-cover, c. These are used to 
determine the surface heat flux, Q*, which is a sum of the four components 
of the flux: 

Q* = Qsw(t,<f>,c) + Qlw(T*,T$,c) +Q E (e s (T a ),e s (T ),v) +Q H (T a ,T ,v) 

(9) 
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which are, from left to right, shortwave, longwave, latent and sensible. The 
functi onal dependencie s for ty pical parameterizations are given (see, for in- 



stance 



Yu and Lindsavl (|2003l )); To is the surface temperature and e s is the 
saturation vapour pressure. The shortwave flux is calculated primarily from 
geometrical considerations. The local solar zenith angle is calculated as 
follows: 

_ / cos 6 cos 6 cos A — sin 6 cos 5 sin A \ 
p = cos 1 = (10) 

\ v cos 2 5 sin 2 A + cos 2 A / 

where phi is the latitude, A is the Earth's tilt, 9 is the time of day as an angle 
(angle of rotation of the Earth around its axis) with 6 = representing noon 
while theta = ir is midnight, 5 is the day of the year (angle of revolution 
about the sun) with 5 = being the winter solstice. 

To return the flux, we multiply the cosine of the local solar zenith angle 
with the solar "constant," S, the cloud cover and the albedo, a: 

Qsw = Scos/3(l - 0.62c)(l - a) (11) 

Spatial derivative were calculated using finite-differencing w hile the equa- 



tion was integrated in time with a fourth-order Runge Kutta (IPress et al 



1992j). The simulation was initialized with a linear temperature profile, 



starting at free zing at the botto m 300 m which was used as the fixed bound- 



ary condition (jRist et all 12002) - It was then spun up by repeated forcing 
with four years of National Center for Environmental Prediction (NCEP) 
reanalysis data. 

2.3 Atmospheric model 

To correct for atmospheric i nfluenc e s we use the parametrised model from 



Wentz and Meissnerl (j2000l ); IWentzl (jl997l ). In the Wentz-Meissner model, 



atmospheric transmissivity and atmospheric components of upwelling and 
downwelling radiation are derived from column water-vapour and cloud wa- 
ter path based on fitted polynomials. The upwelling and downwelling bright- 
ness temperatures are modelled as equivalent temperatures assuming an 
isothermal atmosphere. 

The total downwelling radiation from the atmosphere will be given as 
follows: 

T ba 4= (1 - T a )T d + TT bsky (12) 

where r a is the atmospheric transmissivity, is the equivalent atmospheric 
temperature for the downwelling case and is the cosmic background 
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Figure 1: AMSR-E brightness temperatures at 6 GHz over Lake Superior 
during the winter of 2003. Dots are measured, circles are model results, both 
uncorrected (open) and corrected for atmospheric influences (open circles). 



radiation. The upwelling radiation will be given as: 

T ba t= r a T x t +(1 - t)T u (13) 

where T% f is the modelled surface brightness temperature from the solution 
of Equations (P) and (|2|) and T u is the equivalent atmospheric temperature 
for the upwelling case. If the upwelling and downwelling temperatures are 
equal, = T u , (they are typically very close) then we can model the at- 
mosphere by simply adding an extra layer (the atmosphere) to the ice RT 
model ([1]) and ([2]) and setting the topmost reflection coefficient (interface 
between the atmosphere and space) to zero. 

3 Results 

3.1 Simulation of lake ice 

AMSR-E measurements were collected over the middle of Lake Superior (in a 
radius of 20 km from the coordinate 87 degrees E longitude and 48 degrees N 
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Ta ble 1: Table of atmospheri c correction coefficients calculated according 
to Wentz and Meissner ([2OOO). See Equations (fT2j) and (fT3l) for their use. 
Column water vapour was set at 20 kg/m 2 and cloud water path at 0.1 mm. 
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Figure 2: Measured AMSR-E brightness temperatures at 10 GHz over Lake 
Superior during the winter of 2003. A clustering algorithm has been used 
to separate lake ice from open water points, with the triangles being ice and 
the diamonds open water. 
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Figure 3: AMSR-E brightness temperatures over Lake Superior during the 
winter of 2003. Lines show model results while points with error bars are 
measured. Top is ice, bottom is open water. 
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latitude) for the winter of 2003 during which there was much lake ice. These 
are compared with brightness temperature simulations for the 6 GHz case in 
Figure [TJ A constant ice temperature of zero degrees Celsius was assumed. 
The water temperature was likewise assumed constant at fre ezing and its 



comp lex permittiviy modelled with a Debye relaxation curve (jUlabv et al 
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Measured and modelled values are compared in a brightness temperature 
space consisting of the vertically polarised brightness temperature versus 
the polarisation difference. The cluster of dots on the bottom right is open 
water while the cluster on the top left is lake ice while the larger circles 
are the model results, both corrected for atmospheric effects and uncor- 
re cted. Correction coefficients are shown in Table [TJ calculated according 
to Wentz and Meissner ( 2000l ) assuming total column water vapour of 20 
kg/m 2 and an ice water path of 0.1 mm. 

Simulations are performed for three different ice thicknesses. Fresh water 
ice is almost a perfect dielectric - because of the low value of the imaginary 
part of the permittivity, ice thickness will have less of an effect on the final, 
emitted brightness temperature. Since it is lake ice, we expect it to be 
relatively thin, certainly less than 2 m. 



We use a clustering algorithm ba s ed on kernel density estimation (IMichie et al 



19941 : iTerrell and Scottl . Il992l : iMillsl . [2009) to group the open water and 



ice 



points and separate the two classes from one another as demonstrated in 
Figure [2j Averages for each cluster, along with tolerances calculated from 
the standard deviations are plotted as a function of frequency in Figure [3] 
along with model results. 



3.2 Simulation of Antarctic icepack 

Icepack simulations are a bit more involved because the density of the ice 
varies with depth and because we have to model the internal temperature 
using the thermodynamice model described in Section 12.21 5000 measure- 
ments were chosen at random from AMSR-E swath data in a broad area 
between and 120 degrees East longitude 75 and 85 degrees South latitude 
for the year of 2007. A thermodynamic model forced by four years of NCEP 
data was run for each of the points and then fed to the radiative transfer 
model. A sample temperature profile is shown in Figure HI Results are 
shown in Figures [5] through 



9 



220 230 240 250 260 270 280 
T [K] 



Figure 4: Antarctic ice sheet temperature profile for 27 November 2007, 
14:42:41 UT at 102.97° E, 76.60° S. Profile has been modelled from four 
years of NCEP reanalysis data using a thermodynamic model. 
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Figure 5: Radiative transfer simulations of AMSR-E brightness tempera- 
tures over Antarctic ice pack for the 6 GHz channel. 
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a. Antarctic ice sheet Tbs, 10GHz 
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Figure 6: Radiative transfer simulations of AMSR-E brightness tempera- 
tures over Antarctic ice pack for the 10 GHz channel. 
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a. Antarctic ice sheet Tbs, 18GHz 
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Figure 7: Radiative transfer simulations of AMSR-E brightness tempera- 
tures over Antarctic ice pack for the 18 GHz channel. 
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4 Discussion and conclusions 



The plane-parallel radiative transfer model is able to simulate radiance mea- 
surements from the AMSR-E satellite instrument well, although not per- 
fectly. For the case of lake ice, brightness temperatures are over-estimated 
for the 6 GHz channels, slightly under-estimated for the 10 GHz vertical 
channel slightly over-estimated for the 10 GHz horizontal channel and badly 
under-estimated for the 18 GHz channels, especially in the vertical polariza- 
tion. For higher-frequencies channels, discrepencies are easily accounted for 
by the lack of scattering in the model. As you move up the microwave spec- 
trum from low to high frequency, simple absorption and reflection processes 
become less and less important, while scattering begins to dominate. The 
tendency of scattering, especially within the ice sheet (as opposed to surface 
scattering) is to lower the emissivity along with a complementary increase in 
the polarization difference. Indeed, the measured signal shows an increase 
in the polarization difference which is absent from the modelled signal. In 
fact, the polarization difference over the ice is larger for the measured signal 
over all frequencies. This may be accounted for in part by the presence of 
clouds which were not screened for when selecting measurement pixels. 

Similar discrepencies show up in the open-water signal. In addition to 
scattering caused by clouds and rough surface conditions, the difference may 
also be caused by water temperatures above freezing as the temperature 
of the lake-water was not varied. Sadly, these are still not point-by-point 
validations and with only six data poi nts provide o nly w eak validation of 
the plane-parallel RT model. Recently, Kang et al. ( 2010l ) have derived ice 
thickness from AMSR-E measurements over the Great Bear and Great Slave 
Lakes. This is a statistical model, however, so it would be good to confirm 
the relationship based on more physical considerations, hopefully paving 
the way for more physically-based ice retrieval methods over the ocean, 
just as exist for retrieval of atmospheric parameters. Preliminary results 
suggest that, except for melting ice, the RT model does an excellent job of 
predicting AMSR-E radiances based strictly on thickness. As in the current 
work, ice tempera ture was held constant. Since the lake ice thicknesses in 
Kang et alJ (l2010h were generated using a thermodynamic ice growth model 



temperatures within the ice sheet were returned as a by-product. The author 
hopes to publish a more complete study using these to drive the emissivity 
model. 

For a truly point-by-point analysis, we turn to the model of Antarctic 
icepack. Here the magnitudes are far off, with the model severely over- 
estimating the actual values, however the correlation, while not excellent, 
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Table 2: Biases of the radiative transfer emissivity model compared with 
AMSR-E data over Antarctic icepack. Linear-regression slope and offset 
parameters are presented for all six channels. 



Channel 


Slope 


Bias 


6v 


0.180 


-39.2 


6h 


0.106 


-62.2 


lOv 


0.157 


-40.7 


lOh 


0.086 


-63.3 


18v 


0.239 


-42.3 


18h 


0.148 


-61.4 



is nonetheless strong and significant. As before, the bias is easily explained 
by scattering as gla cial icepack tend s to be very granular with grain sizes 



as large as 15 mm (jRist et al.1 . [2002]). Thus scattering will strongly affect 



even the lowest, 6 GHz channel. Biases are shown in Table [2] which lists 
regression slope and offset values for all six channels. As expected, biases 
are different for the vertical versus horizontal channels in keeping with the 
scattering explanation. Considering that the only input parameters varied 
in the emissivity model were the temperatures, and these were derived in- 
directly from a thermal conductivity model driven by reanalysis data, the 
results seem quite good. 

It is instructive to compare the RT results with a purely statistical model. 
Therefore, 1000 temperature profiles were selected for use as "training data" 
and correlated with the AMSR-E brightness temperatures using the top ten 
(10) singular vectors. Predicted brightness temperatures for the other 4000 
measuremnt pixels are presented in Figures l8l through 1101 along with AMSR- 
E brightness temperatures for comparison. Although the large biases are 
absent, the results are remarkable both for the similarities in accuracy as 
well as in the shape of the scatterplots. This is hardly surprising since the 
RT model is practically linear, with the only non-linearities arising from the 
weak dependence of real permittivity on temperature. 

The weighting coefficients for the RT versus statistical model are com- 
pared in Figure PTlT a) and lllf b). For the statistical model, simulated bright- 
ness temperatures are simply a linear combination of the temperatures of 
the icepack. To account for the weak nonlinearity in the RT model, weights 
are calculated with a numerical derivative, averaged over 500 trials: 
standard deviations are also shown in the figures. 
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Figure 8: Statistical model of AMSR-E brightness temperatures over 
Antarctic ice pack for the 6 GHz channel. 
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a. Antarctic ice sheet Tbs, 10GHz 
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Figure 9: Statistical model of AMSR-E brightness temperatures over 
Antarctic ice pack for the 10 GHz channel. 
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a. Antarctic ice sheet Tbs, 18GHz 
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Figure 10: Statistical model of AMSR-E brightness temperatures over 
Antarctic ice pack for the 18 GHz channel. 
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q. Weighting functions h pol. 
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Figure 11: Temperature weighting functions for Antarctic ice sheet for hor- 
izontal (a) and vertical polarizations (b). The weighting functions for the 
radiative transfer simulations represent an average over fifty (50) AMSR-E 
measurement pixels randomly selected from the test set. The averages are 
shown by the solid line while the dotted lines enclose the standard devia- 
tions. The broken lines show the weighting coefficients for the statistical 
model. 
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The two types of weighting functions are quite similar, although the 
statistical weighting functions are shifted considerably upwards. They also 
allow for negative weight. Thus, scattering means that the ice pack is ef- 
fectively much more opaque than our non-scattering model suggests. This 
provides a further explanation for the too-high brightness temperatures re- 
turned by the RT model since ice temperatures tend to be higher further 
down in the icepack-see Figured! Nonetheless, the penetration depth of the 
microwaves is quite deep: as much as 50 m for the lowest frequency. This 
means that radiometer measurements of glacial icepack return information 
about temperatures within the ice sheet, which are reflective of historical 
temperature. 

The interaction of electro-magnetic radiation with matter is a very fun- 
damental problem. Most operational ice retrieval algorithms are based pri- 
marily on empirical parameterizions. While the model under discussion 
is somewhat ad- hoc, nonetheless it is based on real physics. Such physi- 
cal models are invaluable for ice retrieval efforts, both for use as forward 
models and to retrieve physical characteristics of ice sheets. Most likely 
parameters to be retrieved include ice temperature, ice thickness and in 
particular, complex permittivities within the ice sheet which are reflective 
of the ice composition and micro-structural properties. Retrieving fresh- 
water ice thickness appears to be quite easy, however retrieving salt-wa ter 
ice thicknesss is proving to be quite difficult ((Mills and Hevgsterl . fioilbl lah . 

Ice, particularly saline ice, is a highly complex composite. Because of 
this complexity, determining its electromagnetic properties -how it interacts 
with radiation-is very difficult. More work needs to be done to understand 
how ice emits, absorbs and scatters electromagnetic radiation. 
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